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A phase space theory for population dynamics in Ecology is presented. This theory applies for a 
certain class of dynamical systems, that will be called Al-systems, for which a conserved quantity, the 
Al-function, can be defined in phase space. This A^-function is the generator of time displacements 
and contains all the dynamical information of the system. In this sense the A^-function plays 
the role of the hamiltonian function for mechanical systems. In analogy with Hamilton theory we 
derive equations of motion as derivatives over the resource function in phase space. A Al-bracket 
is defined which allows one to perform a geometrical approach in analogy to Poisson bracket of 
hamiltonian systems. We show that the equations of motion can be derived from a variational 
principle over a functional J of the trajectories. This functional plays for Al-systems the same 
role than the action S for hamiltonian systems. Finally, three important systems in population 
dynamics, namely, Lotka-Volterra, self-feeding and logistic evolution, are shown to be Al-systems. 

PACS numbers: 87.23. Cc 87.10.-|-e 87.23.-n 



Conservation of energy modeled by differential geom- 
etry is the cornerstone of the phase space formulation of 
classical mechanics p^ . There exist a wide machinery of 
analytical methods developed during the past two cen- 
turies for such systems |5]. For instance, the Hamilton 
formalism resorts on a Hamiltonian function defined on 
phase space which contains all the dynamical informa- 
tion of the system. Once specified, we can generate upon 
differentiation equations of motion for the system. 

The situation is much more precarious in population 
dynamics in Ecology. Ecological systems are overwhelm- 
ingly complex systems [3]. However, in some cases they 
can be described by simplified phenomenological rate 
equations that describes the net effect of the environ- 
ment and the interaction between densities of species [3] . 
However, a systematic approach allowing one to derive 
equations of motion from general principles is still lack- 
ing. The present paper is a step towards this goal. In 
fact, in some cases a mutual resource function is con- 
served and keep constant during the evolution of the in- 
teracting species. The question is then, can be develop 
a phase space formalism that expresses the conservation 
of this quantity? 

We show in this paper that there exist a wide class 
of systems in ecology, that will be call A^-systems, for 
which population dynamics can be described by a phase 
space theory analogous to the hamiltonian formulation of 
classical mechanics. We exploit this analogy in order to 
use the elaborated analytical methods developed for the 
latter. Hamilton function, Poisson bracket, variational 
principles, action functional and Lagrangian, all of them 
useful tools for hamiltonian systems will find their coun- 
terparts in the theory. 

Finally, we show that three important system in pop- 
ulation ecology are A^-systems. They are the Lotka- 
Volterra predator-prey system, the self-feeding interac- 
tion and the logistic evolution. These systems are also ap- 



plied in many other disciplines outside ecosystems, which 
sets forth the wide range of applications of the A/(-system 
formalism. 

We start with the definition of A/J-systems. A J\4- 
system is a dynamical system defined by a given resource 
function A4{q,p ; i), and equations of motions given by 
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Not every conceivable motion in phase space F = iq,p) 
is a A4-dynamical system. In order to obtain a necessary 
condition let us consider the generic dynamical system 

q = gi{q,p), 

P'^g2{q,p), (2) 

where gi,g2 G -7^(r) are two arbitrary functions in phase 
space. The existence of a continuous A4-function requires 

d'^M _ d^M 
dqdp dpdq 
Eq. ^, together with and Q yields 

9]_ _^ g2 ^ dg^ ^ dg2 ,^. 

q p dq dp ^ 

which is a necessary condition for the dynamical system 
(pi) to be a A^-system. If this is the case we say that 
51,52 are the components of a A^-vector field. 

According to Eq. (fTl), the resource function Ai char- 
acterizes the dynamics of the system. The A^-systems 
are not hamiltonian systems, since they do not obey the 
hamiltonian canonical equations of motion. However, 
they possess a close relationship with each other. In fact, 
consider the change of variables 



q' = ln(q/cr), 

p' = Hp/0, 



(5) 
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where a, ^ are two positive constants. Using Eqs. (fTl) it 
is easy to check that they obey the hamihon canonical 
equations [1] 



(6) 



where H{q',p') ~ A4{q{q'),p{p')) is the hamiltonian 
function. Thus, Eqs. (Is]) aUows one to construct an as- 
sociated A^-system starting from a generic hamihonian 
system. Conversely, note that these relations comprise 
the Hamilton phase space A = M^ into the first octant 
M."^ = (q,p e R^q^p > 0). Thus, the class of A^-system 
on r = M^ has an associated hamiltonian system on A. 
We will discuss later the importance of this system sub- 
class in ecology and game theory. 

For many dynamical systems Ji gives their mechanical 
energy E. The resource function A4 play for A^-systems 
the same role than the function Ji for hamiltonian sys- 
tems. In order to see this, we introduce a geometrical 
approach analogous to the Poisson bracket of hamilto- 
nian system. The Poisson bracket (PB) is defined as [1 



{f,9}PB^{d,f)n^'{dkg) 



(7) 



where f,g(z J^{T) are two arbitrary functions on phase 
space, i7 is tlie symplectic matrix 



n 



1 
-1 



(8) 



d 



and we have defined the derivative operators di 
where r/i = q, 1^2 = P- We follow the Einstein rule of 
summation over repeated indexes. In terms of the PB, 
Eqs. (|6| can be written as 



iji = {f?i,'H}pB. 



(9) 



Moreover, the evolution of any dynamical variable 
/(q',P;0 G r is given by 



f^{/,H}pe 



5/ 
df 



(10) 



According to (10 1, H is the generator of time displace- 
ments of the system. Due to the antisymmetry of the PB, 
the time evolution of the hamiltonian is H = ^ . This 
reflects the conservation of energy for time-independent 
hamiltonian systems. 

We now define the Af-bracket (MB) as 



[f,g]MB = id,f)n^Hdkg), 



(11) 



where we have defined the derivative operators di = 
otA — r. It can be shown that MB (as the PB) is bilinear, 
antisymmetric and satisfies the Jacobi identity. These 



are the defining properties of a Lie algebra [S] . The func- 
tion space J'(r) is a Lie algebra under the action of the 
MB. MB also satisfies the product rule 



[/,5^]mB = [/,ff]MB/i + 5[/i^]MB- 



(12) 



Thus, MB is a kind of derivative operator [6]. As a matter 
of fact, Eqs. (IT]) can be recast in terms of the MB as 



m = V\i,M\MS,- 



(13) 



The function At G -^(L) is the generator of time dis- 
placements on the system, since the time derivative of 
any dynamical variable / along the motion is 



dt 



[/,-M] 



MB 



df 



(14) 



dM 
dt 



In particular, choosing Ai — f we obtain A4 
the resource function of the system is conserved if it does 
not depend explicitly on time. The equation A4(q,p) — R 
for Al-systems plays the same role than the equation 
'H{q,p) = E for hamiltonian systems. 

Hamiltonian systems admit a variational formulation. 
In fact, dynamical systems in classical mechanics mini- 
mize the action functional [T] 



L{q'{t),q'{t),t)dt (15) 



where L is the Lagrangian function of the system 

L^p'q'-H (16) 

It can be shown [1 that Hamilton canonical equations 
given in (p| can be derived from the variational principle 
dS = 0. With the help of Eqs. ([5|, we can now construct 
a variational principle for the associated Al-system. To 
this end we define the A^-Lagrangian 



C=^\nip/0-M, 

q 



and the A^- Act ion 



J^ / Cdt. 

Jta 



(17) 



(18) 



It is easy to check [7] that the imposition of the extreme 
condition to the Al-Action, i.e., 5 J = 0, renders the 
canonical equations of motion for Af-systems given in 
Eqs. (fTl). Thus, among all the possible trajectories that a 
A(-system could undertake, the actual motion {q{t),p{t)) 
is the one that makes Eq. (jTsl) minimal, i.e., SJ^ — 0, 
where 6J^ is the variation of J7 over different trajectories 
connecting {q[to),p{to)) and {q{ti),p{ti)). 

J and £ play for the Al-systems the role that the me- 
chanical action S and the Lagrangian L play for hamil- 
tonian systems. The definition of H retains a certain 



amount of ambiguity |S]. As a matter of fact, we can 
construct H in the form 



£ = qa +P/3 - 7W, 



(19) 



where a(q,p), P{q,p) are generic functions in phase space 
(dynamical variables) satisfying 



dq 



1 
qp 



da 
dp' 



(20) 



The main goal of population dynamics is to describe 
the evolution of the population of different species in- 
teracting in an ecological system. We show now that 
Al-systems comprises a wide class of important system 
in ecology. Let us start with the Lotka-Volterra (LV) 
model [H [To] describing an ecological predator-prey (or 
parasite-host) system. The LV two-species model is given 
by the equations 

"^^ (ai - b,V) U, 



dt 
dV 
lit 



{~a2 + b2U)V, 



(21) 



where ai is the the growth rate of prey C/, 6i the rate at 
which predators V destroy the prey, 02 the the death rate 
of predators and 62 the rate at which predators increase 
by consuming prey. We take the adimensional variables 
proposed by Hsu [TT] 



0.2 



U, p- 



h 

0.2 



y. 



a2t, C 



02'' 



(22) 



to reduce the Lotka-Volterra (LV) system to 
dq 

dp 
dr 
In spite of non-linearity of the equations of motion 



(e-p)9, 



(23) 



given in (23 1 an analytical solution is possible, namely. 



q p^ ^ e' 



where R is an integration constant, 
define 



(24) 



This takes us to 



A^ ~ (^Inp — p + hiq — q, 



(25) 



since A^ = i? is a constant of motion. Moreover, it can 



be rapidly checked that the insertion of (25) into Eqs. 
^ yields the LV equations given in ([23]) [12 . Thus LV 
is a A^ dynamical system. 

This ecological interpretation of the Al-function as 
a mutual abstract resource shared by the interacting 
species is enforced by considering a simple self-feeding 
system, modeled by the equations of motion 



q 

P 



qp, 

-qp. 



(26) 



consider direct disease transmission [T3] in a total pop- 
ulation oi q + p individuals where q are the number of 
infected individuals and p the number of individuals who 
are susceptible of infection. In this simple system there 
is an obvious conserved quantity: the total population 
given by the initial condition R = qo + Po- R can be 
interpreted as a mutual global resource being distributed 
between the "interacting species" q and p that remains 
conserved during the process. The corresponding M- 
function is then 



M = q + p. 



(27) 



In fact, it can be checked rapidly that the equations of 
motion (26) are obtained by inserting (27) into Eqs. (II]). 
Consider now the A^-function for a two-species 
Malthusian system 



M 



-"flnq + fclnp, 



(28) 



SO q 



hf 

qoc''^ p 



PoG^ ■ Consider now the coupling 



Ai = —7 Inq + klnp ' 



1 P 
win—, 

Ps 



(29) 



where Ps > Po is a saturation constant. The canonical 
equations of motion given in Eqs. ([I]) equipped with ( 29 ) 
yield 



p = '-tp + fiqpln{p/ps). 



(30) 
(31) 



Eq. ( 30 ) defines logistic evolution, where k is the grow- 
ing rate for the population q. The mortality rate ^q is 
density dependent, making the logistic evolution a sim- 
ple model to describe the phenomena of saturation of 
population growth in time [ 151 I16j . 



Now we turn to the analysis of Eq. (31 ). Remarkably, 
Eqs. (30) and (31) have both analytical solutions, p is 



the amount of food supply feeding the system, and thus 
being depleted by q. We take 7 = and the change of 
variables t — kt, p — p/ps and q ~ q/qs where qs is the 
saturation value qs — k/ji. Eqs. (30)-(31) reads now 



dq 

dr 
dp 

dr 



q{i~q), 

qp In p. 



(32) 
(33) 



An analytical solution of Eqs. (32)-(33) can be found 
in the form 



(34) 
1 + - I , (35) 



e^ -1 
p = exp 



qo 



where a < is an integration constant fixing po = e"' ''° 
Therefore, a depends on the distance between po and Ps 




FIG. 1. Logistic solutions for q (thin solid line) and p (dashed 
line) as a function oi t — kt for a = —0.03, Qo ~ 0.14, g^ = 
1.41, Ps = 1.55 and k — 10. The resource function (thick solid 
line) remains constant in the evolution A4 — k{\nps — D). 



It is easy to check that the solutions obtained make Ai 
constant in time. 



The solution (35 1 has two stationary limits, p = Ps and 



the asymptotic limit p = 0, since both limits make the 
right hand side of Eq 
and (a < 0, T - 



(33) vanish. They correspond to 
respectively. The parameter 



a = U ana (a < U, r — >■ cx)) 
a measures how fast the food supply is depleted; giving 
the width of the decay as A = ln(l — 1/a). This is 
shown in Fig. IT] where Eqs. (34) and (35) are plotted for 
a = -0.03. 

An explicit form of p as a function of q can be obtained 
after elimination of the time variable between Eqs. (34) 



and ( 35 1 , namely 



P 



-D/(l~q) 



(36) 



where D = a{l — l/qo)- The above equation provide us 
of the phase space of the system. 

The analogy between A^-systems and hamiltonian sys- 
tems can be exploited further. In this paper, a geomet- 
rical approach has been developed based in the intro- 
duction of a A^-bracket which, as the Poisson bracket of 
hamiltonian systems, has the structure of a Lie algebra. 
The 7W-bracket establishes a bridge than can be used 
to generalize the breadth of analytical methods in phase 
space developed for the Poisson bracket [2] (Noether the- 
orem connecting symmetry and conservation laws, Liou- 
ville theorem establishing conservation in time of phase 
space volume for ensembles, etc). 



The author thanks ecologists Antonio Perez and F. 
Garcia Novo for pointing out some helpful bibliography. 
J.M-L is supported by a Return Program from the Con- 
sejeri'a de Educacion y Ciencia de la Junta de Andalucia 
in Spain. 



[1] 
[2] 



[4] 
[5] 
[6] 



[7] 



H. Goldstein, C.P. Poole and J.L. Safko Classical Me- 
chanics 3 edition (Addison Wesley, NY, 2002). 
J.V. Jose and E.J. Saletan, Classical Dynamics, Cam- 
bridge University Press (1998). 

[3] M. Begon, J.L. Harper and C.R. Townsend, Ecology 
(Blackwell Science, Cambridge, 1996). 
Frank C. Hoppensteadt, Mathematical Methods m Pop- 
ulation Biology (Cambridge University Press, 1982). 
D. H. Sattinger y O. L. Weaver, Lie Groups and Algebras, 
(Springer- Verlag, New York, 1986). 

Another well known example of C* Lie algebra is the 
one associated to the Moyal bracket, generating time dis- 
placements of the Wigner function. This geometrical ap- 
proach yields the phase space formalism of Quantum Me- 
chanics given in J. Moyal, Proc. Camb. Phil. Soc. 45, 99 
(1949). 

For instance start form the modified Hamilton principle 
in phase space, which imposes Euler-Lagrange equations 
to £, as can be seen in [2]. 

[8] This ambiguity is inherent to the Hamilton modified prin- 
ciple. In fact, it also applies to the definition of L, as can 
be seen in [l]. 

Lotka, A. J. 1925. Elements of physical biology. Balti- 
more: Williams and Wilkins Co. Reissued as Elements 
of Mathematical Biology. Dover, New York, 1956. pp 88- 
90. 

Volterra, V. 1926. Mem. R. Accad. Naz. dei Lincei. Ser. 
VI, vol. 2. A complete translation can be found in Ap- 
plicable Mathematics of Non-physical Phenomena, F. 
Oliveira-Pinto and B.W. Conolly. John Wiley and Sons. 
1982. pp 23-115. 

S.-B. Hsu, J. Math. Anal. Appl. 95, 428-436 (1983). 
[12] This fact was already found in the old papers [T^ for the 
particular LV system. However, our algebraical approach 
allows us to construct a phase space formalism for general 
Al-systems. 

E.H. Kerner, BuU. Math. Biophys., 19, 121-146 (1957), 
E.H. Kerner, BuU. Math. Biophys., 21, 217-255 (1959), 
E.G. Leigh Jr., Zoology, 53, 777-783 (1965). 
J.H. Vandermeer and D. E. Goldberg Population Biology. 
(Princeton University Press, New Jersey, 2003). 
C.J. Krebs, Ecology: the experimental analysis of distri- 
bution of abundance. Harper and Row (1985). 
For a list of assumptions implicit in the LV and logis- 
tic equations see for instance R.W. Poole, Quantitative 
Ecology, (McGraw-Hill Series in Population Ecology, NY, 
1974), pgs 148 and 63 respectively. 



[9] 



[10] 



[11] 



[13] 

[14] 
[15] 
[16] 



